Typeset with jpsj2.cls <ver.l.2> 



Short Note 



Calculation of the Self-energy of Open Quantum Systems 

Keita Sasada 1 * and Naomichi Hatano 2 ^ 

1 Department of Physics, University of Tokyo, Komaba, Meguro, Tokyo 153-8505 
2 Institute of Industrial Science, University of Tokyo, Komaba, Meguro, Tokyo 153-8505 

KEYWORDS: Open quantum system, Self-energy, Resonant state 

The electronic conduction in mesoscopic systems has been studied extensively in recent 
years. A theoretically interesting feature of the problem is the fact that the system in question 
is an open quantum system with semi-infinite leads. The open quantum system intrinsically 
has resonant states, which can strongly affect the electronic conduction. 1 

A popular way of treating the semi-infinite leads is to contract the leads to the self-energy. 
The self-energy of leads is a useful way of computing the conductance as well as obtaining 
resonant states. In this note, we propose a new method of calculating the self-energy of the 
leads. The self-energy T,(E) was originally defined in 2 

{A E-H + l 5 W) = ^ E-{H c + nE)) W) (1) 
for sites x and x' inside the central conductor, where H c is the Hamiltonian of the central 
conductor and H is the total Hamiltonian including semi-infinite leads attached to the con- 
ductor. The self-energy has been calculated by various methods. The method that we present 
here is much easier than previous methods. The main claim of this note is that the self-energy 
is equivalent to the boundary conditions for resonant states. 

We consider the Hamiltonian of a conductor with semi- infinite leads attached to it: H = 
H c +Y^ a H a , where H c is a one-body Hamiltonian of a finite-size conductor, while H a describes 
a semi-infinite lead given by the tight-binding model 

oo 

H a = -tY^ (\x a + l}(Xa| + \x a )(x a + 1|). (2) 
x a =0 

This includes the hopping between a site x a = on the conductor and the lead a. (Note that, 
if we have hopping between the conductor and a lead with the amplitude different from —t, 
we include it in H c .) 

Equation (1) suggests that the eigenvalues of the effective Hamiltonian H c r{E) = H c + 
are the poles (bound states and resonant states) of the total Hamiltonian H on the 
complex E plane. Therefore, we seek discrete and generally complex eigenvalues E n of resonant 
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states and bound states of the whole system: 

H\4> n ) = E n \4> n ) and $ n \H = E n $ n \. (3) 

The eigenfunctions are bi-orthogonal: (tp n \'4'm) = $nm- The eigenvalues E n are related to the 
corresponding eigen- wave-number k n , which is also generally complex, through the dispersion 
relation E n = —2tcosk n . The eigen-wave-number k n is on the upper-half plane for the bound 
states and on the lower-half plane for the resonant states. 

It is known that the resonant states as well as the bound states can be found by requiring 
the boundary conditions (x a \i/) n ) oc e thnXa for x a > in the leads. 3 In other words, the discrete 
states satisfy the boundary conditions 

(x a + l\ip n ) = e ikn (x a \i/> n ) for x a > 0, (4) 

where $tk n > 0. The boundary conditions (4) transform the Schrodinger equation 

(x a = 0\H c \lp n ) - t(x a = l\ip n ) = E n (x a = 0\ip n ) (5) 

to 

(x a = 0\H c \^ n ) + V^\E n )(x a = 0|V„) 

= E n (x a = 0|^n), (6) 

where 

V £\E) = -te ik (7) 

is the energy-dependent effective potential. 

We claim that the self-energy of the lead a is nothing but the effective potential: 

£ (Q) (£) = v£\E)\x a = 0)(x a = 0|. (8) 
The total self-energy is the sum over the leads: = ^^(E). The effective potential 

a 

is rewritten in terms of E as 

by using the dispersion relation E = —2t cos k. Note that we choose the branch $>V^ < for 
the retarded Green function. Equation (9) is indeed equivalent to the expression obtained by 
other methods. 2 

Let us now demonstrate that the present method is easily generalized to other types of 
leads such as N-leg ladder and carbon nanotube. Hereafter, we drop the lead index a for 
simplicity. First, we calculate the self-energy of a lead of TV-leg ladder (Fig.l): 

oo N 

^ladder = ~ t ^ J ^2 (\ x + 1 ,V}( x ,y\ 
x=0 y=l 

+\x,y + l){x,y\+c.c). (10) 
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Fig. 1. Leads of the form of ladders are attached to the central conductor. 



We first diagonalize ladder m the V direction and obtain the conduction channels 
{4>j{y) \j = 1,2,-- - ,N}, where 



>j(y) = sin 



N + l, 



N 
\ V"=l 



sin 



2 Jny" 



N + l 



(11) 



Each channel has the dispersion relation E = — 2icos kj + ujj, where Uj = —2tcos(2irj/N). 
Each channel yields its effective potential of the form Eq. (7), or 



V®{E) = -te*> = - - 



(12) 



(13) 



The self-energy of TV-leg ladder is given in the N x N matrix form 

N 

(SjadderC^))^ = ^ ^)Kff 0/)*- 

i=i 

The result is equivalent to the one obtained in Ref. 4. 

Second, we calculate the self-energy of a lead of (n,0) zigzag carbon nanotube attached 
to the conductor as in Fig. 2, where n is the chiral number. The Schrodinger equation of the 
zigzag carbon nanotube flzigzag I^a/b^O) = -^I^a/b(%)) ytelds the dispersion relation of the 
jth channel as 



E= ±t \hi 



±M/l±4cos 



V3kj 



— - cos h 4 cos z — , 

2 n n 



with 



h k = eVs +2cos— e ^ 

3 n 



(14) 



(15) 



where the first Brillouen zone is \kj\ < and its wavefunction on the A and B sub-lattices 

as 

h* 

k -u 

\h kj \ 9 ° yUh (16) 
(x,y\^(k j ))= e^-(y), 



{x,y\ip^(kj)) 
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where 4>j(y) = e l ~^~ y i 
of the whole system 



where 4>j(y) = e l " jy /\/n. The boundary conditions (16) transform the Schrodinger equation 



(x = 0,y\H c \^(kj)) - t(x = 1/V3,y\^(kj)) 

= E{x = 0,y\l£(k j )) (17) 

to 

(x = 0,y\H c \^(kj)) + v}ffJ g (E){x = 0,y\^(kj)) 

=E(x = 0,y\i>±(k j )), (18) 
where the effective potential of the jth channel is given by 

V<j™W - (19) 



E 2 + t 2 — X 2 ± iyj (2tX j ) 2 - (& -t 2 - A 2 ) '' 



with \j = 2tcos7rj /n. Hence we obtain the self-energy of an (n,0) carbon nanotube in the 
n x n matrix form 

n 

(EWtf) W = E MVB)V$&{E)<l>Mr> (21) 

J'=l 

where j/a and ye are coordinates on the A and B sub-lattices, respectively, which are indicated 
in Fig. 2. The result (21) is indeed equivalent to the one obtained in Ref. 6. 

When the A sub-lattice, instead of the B sub-lattice, is in contact with the conductor, we 
obtain the self-energy in the form 



V ^zigzag \^ JJy A) y' A - / sVjWAJV 

with 



PWtf) W = E MvaK^E)^)* (22) 



V ir ' A) (E) 

zigzag V / 



E 2 - t 2 + X 2 ± iJ (2tA,) 2 - (i? 2 - 1 2 - A 2 ) 

— . (23) 



2£ 
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Fig. 2. A lead of the zigzag carbon nanotube. The upper and lower edges satisfy the periodic boundary 
conditions. 
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